DO-TH 95/15 
RAL-TR-95-043 
August 1995 



Spin-dependent non-singlet structure functions 
in next-to-leading order 

M. Stratmann, A. Weber* 

Universitiit Dortmund, Institut fiir Physik, 
D-44221 Dortmund, Germany 

W. Vogelsang 

Rutherford Appleton Laboratory 
Chilton Didcot, Oxon 0X11 OQX, England 

Abstract 

We study in detail the fiavor-non-singlet component of polarized structure functions in 
the framework of a consistent and complete next-to-leading order {0{as)) analysis. In 
this context, we discuss some important features of the calculation of the next-to-leading 
order corrections. Particular emphasis is put on the Q^-evolution of sum-rules for the 
first moments of the non-singlet structure functions which, as we show, could serve to 
explore SU (2) and SU (3) breaking effects in relations between baryonic /5-decay matrix 
elements and in the proton's polarized sea. Furthermore we make predictions for polarized 
non-singlet structure functions measurable in a conceivable ep collider mode of HERA. 
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1 Introduction 



The spin-dependent structure functions of protons, neutrons, and deuterons have received 
much attention both experimentally and theoretically in the past years. Since the advent 
of the EMC result on the proton's gi{x, {Q'^) = 10.7GeV^), most theoretical studies 
have been focussed on the singlet component of this structure function in order to explain 
its unexpected experimental smallness, hereby assuming that the non-singlet (NS) com- 
ponent is rather well understood. First experimental evidence for this latter assumption 
was provided recently by the confirmation 0, ^ of the Bj0rken sum- rule 0] which relates 
the integrals (first moments) of and g'^. However, this sum-rule, which depends merely 
on the fundamental SU (2) isospin {u ^ d) symmetry between matrix elements of charged 
and neutral axial currents and is therefore expected to hold, does not entirely fix the first 
moment of the NS component, gi^si 9ii since the latter can be written (in leading 
order (LO)) as 

gUs = ^^M^) + YQ^M^) , (1) 

where in terms of the first moments of the polarized (anti) quark densities A 1 {x, Q"^) we 
have 

AAsil) = (^Au + Au- Ad- Ad) {x,Q'^)dx , 

AAsil) = (^Au + Au + Ad + Ad-2{As + As)) {x,Q^)dx . (2) 

The Bj0rken sum- rule is equivalent to 

AA3(1) = F + L) = c/A = 1.2573 ±0.0028 , (3) 

but information on AAg can only be obtained from hyperon /3-decays. Assuming full 
SU{3) symmetry between hyperon decay matrix elements of the flavor-changing weak 
axial currents and the neutral ones, one finds (with F, D taken from ref. |P) 

AA8(1) = 3F-L) = 0.579 ±0.025 . (4) 

This approach has been seriously questioned in ref. |Q, where the suggestion was made 
that SU (3) f symmetry is broken in such a way that only the valence Aq^ = Aq — Aq 
content of Ay43(l), Ay48(l), rather than the full combinations AA^^ll), Ay48(l), enters eqs. 
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(3), (4). In view of these uncertainties and of the fact that the baryonic /^-decays cannot 
tell us anything about gi^^s except for the first moment, it is interesting to examine the 
NS sector of polarized structure functions in more detail in order to find other possible 
experimental clues to A A3, AAg, and the polarized valence densities Am^,, A(i„, hereby 
improving our present understanding of the relation between the first moments of these 
quantities and the F and D values. In this respect, it is necessary to consider not only 
the NS piece of the electromagnetic (e.m.) structure functions (/f " = (^i^'^", but also the 
polarized electroweak structure functions g^, studied in refs. [7-14], which partly are 
pure NS quantities. Since possible measurements of such structure functions are likely 
to be performed at much higher than those relevant for eqs. (3), (4), it is important 
to theoretically understand the Q^-evolution of spin-dependent NS structure functions as 
well as possible. For this purpose, it is necessary to improve the theoretical predictions 
in the NS sector, by performing a complete and consistent next-to-leading order (NLO) 
analysis. All ingredients for this are available as we will see below, and there are some 
features of the NLO corrections which are interesting in themselves. Also, the theoretical 
predictions are much more reliable for the NS sector since it is not plagued by the anomaly 
contribution as is the case for the singlet contributions to polarized structure functions 

The remainder of this paper is organized as follows: In Section 2 we review the main 
LO results on spin-dependent NS structure functions. Section 3 gives a detailed account 
of the determination of the NLO corrections. Section 4 contains the numerical evaluation 
of our results, and in Section 5 we summarize our findings. 



2 Spin-dependent non-singlet structure functions in 
leading order 



The structure functions gi, (73, and g^ appear in the hadronic tensor as (see, e.g., W% ) 

W^. = -ie^upaj^gi + [-g^.u + ^ ) ^3 + ^ K ~ ) K ~ ) ^4 ' 

where we have already replaced — > for the spin vector of a longitudinally polarized 
nucleon with momentum p. In eq. (5), q denotes the momentum of the virtual boson 



probing the hadron. As is well-known, and (74 do not contribute to purely e.m. scatter- 
ing, but appear in (parity-violating) electroweak neutral current (NC) or charged current 
(CC) lepton-nucleon interactions. Therefore, their experimental accessibility may seem 
remote presently. However, they could be measured in 'i/ -scattering off a polarized tar- 
get, and they would certainly play a role in deep-inelastic scattering (DIS) experiments 
at HERA if also the 820 GeV proton beam could be longitudinally polarized All 



relevant cross section formulas for ^up, e^p interactions in terms of gi, g^, and (74 can be 
found, e.g., in refs. |jT2|, |13| and need not be repeated here. As was shown in |12|, the LO 
expressions for the structure functions can be cast into the forms 

9i{n,Q') = WsJAq{n,Q')+Aq{n,Q')) 
^ q 

93{n,Q') = WRA^^i^^Q')-^^i^^Q')) (6) 



9Ain - 1, Q') = J2Rq i^li^^ Q') - ^^(^' Q')) = 2^73(n, Q^) 
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where, as usual, the Mellin-n moments of a Bj0rken-a:-dependent function g{x) are defined 
as g{n) = x'^~^g{x)dx. The coefficients Sq and Rq in eq. (6) depend on the exchanged 
boson, 7*, W^, in the DIS process and can be found in ||l2l. Obviously, for W^- 
exchange (CC interactions), only the quark or the antiquark of a given flavor contributes, 
depending on the charge of the W. To become more specific, we write the various con- 
ceivable structure functions in LO in terms of the NS quark combinations (for eqs. (7-20) 
below we drop the obvious argument (n, Q^) from all quantities) 

Aujj = Au — Au , Ad^ = Ad — Ad , 
AA3 = Au + Au- Ad- Ad , 

AAs = Au + Au + Ad + Ad- 2(As + As) (7) 



and the singlet 



(for / = 3 fiavors) as 



AS = Am + Au + Arf + Arf + As + As 



^r^^-'"-^ = ^AA3 + ^AAg + ^AS (9) 
9?''''' = \[lAu^ + lAd^) (10) 
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for ep NC interactions, where for Qi we have only written the purely electromagnetic 



contribution (which dominates |12] if polarized electrons are used) since the other NC 
contributions do not easily lead to NS quantities, g'^^'^^ has been written only for the 



dominant |12| contribution from '-yZ^ interference^. For CC structure functions ('z/'p 
e^X or e^p — >*z/'x scattering) one has [|l^ 

gl^ = Arf + AM + /i(A)As = ^(Arf,-AM, + AS)-^(l-/i(A))(AS-AA8)(ll) 

gl^ = Au + Ad + MX) As = I {Au, - Ad, + AS) - hi - MX)) (AS - AAs) (12) 

2 o 

g'^" = -Ad + Au- /3(A) As = I {AAs - Au, - Ad,) - ^^ (AS - AAs) (13) 

2 o 

g'/ = -Au + Ad + /3(A)As = {AAs + An, + Ad,) + ^^ (AS - AAs) , (14) 

2 

where we have introduced functions fi{X = Q"^ / {Q"^ + ml)) with /i(l) = 1 (massless limit) 
which take fully into account the effects of the charm mass ttIc in the s c transition. 
In the LO considered in eqs. (11-14), the fi{X) are simply given by the 'slow-rescaling' 
prescription Jl^ which yields xsj = Q'^/'^VI ~^ xsj/X and therefore /i(A) = /3(A) = 
/(A) = A" for the nth moment in eqs. (11-14). The expressions for the rath moment of the 
structure function gi/2x, (74(72 — 1, (5^)/2, are the same |]^, ^, |Tl|, |1^ as the right-hand-sides 
(rhs) of eqs. (10), (13), (14) with, however 0, /4(A) = /(A)/A in the CC case. In this way 
one finds a (slight) violation of the Callan- Gross-like relation (74(71,-1, Q"^) = 2(73(72, Q"^) by 
terms of 0{ml/Q'^) due to the CC s c transitions. Note that eqs. (11-14) can be easily 
seen to receive only corrections of 0{sm'^ 9cml/Q'^) when taking into account the effects 
of Cabibbo-mixing; we can safely neglect these small terms. The structure functions for 
DIS scattering off neutron-targets can be easily obtained by changing the signs of the 
AA3 terms and interchanging Am, Ad, in eqs. (9-14). From eqs. (9-14) we can, e.g., 
construct the following NS combinations 0: 

ep(e.m.) en (e.m.) 1 A /I /i 

^1 -gi = ^A^a (15) 

gr-g^ = Au,-Ad, (16) 

g7~gT = AA3 (i7) 
g^ + gr = -{Au, + Ad,) (18) 



^ As compared to ref. |12 we have dropped a factor (4sin^ Qw cos^ QwiQ^ + M^)/Q^) ^ from the 
normahzation of g'^'^'^ ^ where Qw is the Weinberg angle and Mz the mass. 
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ep (e.m.) , enie.m.) 

+ 91 



9? + 9? 



9 {9I 

2 + /(A) 
/(A) 



9? - 9T 



2 + /(A) 



{9? + 



15/(A) -2 
2 /(A) + 2 



AAg (19) 
(20) 



etc. Besides these relations, the NC g'^'^^ [n, Q"^) in eq. (10) is obviously also an entire 
NS quantity. 

3 Next-to-leading order corrections 



In order to study the evolution of the polarized NS structure functions in NLO it is 
necessary to recall the well-known solution (see, e.g., |1^) of the NS renormalization group 
equation relating the Mellin-n moments of a polarized NS structure function cji^NS = 9ii 
93, 94:/'^x at the input scale Ql and at Q'^ > Q"^: 



9i,Ns{n,Q^' 



1 + 



as{Q^)-as{Q 



An 

2\\ 7^sW/2/3o 



2AClin) + 



2Po 



2Pi 



X 



'«s(<5 



where 



1 + 



2n 



9i,Ns{n,QQ) 
-ACVn)^ Agi,5(^,Q') 



/5i InlnQVA 



MS 



47r /?o In QVA|s Pi {\^ Q2 /a_^^^ 
with the QCD scale parameter A^, /5o = 11 — 2//3, Pi = 102 — 38//3, and 



(21) 
(22) 

(23) 



2Po 



2P'o 



X 



' asm 



(24) 



is the suitable NS combination of polarized quark densities evolved from Qq to via 
the LO (one-loop) and NLO (two-loop) NS anomalous dimensions^] 7^^(72) and Ajl^g^n). 
The precise form of the NLO pieces (AC*, Ajlfg) in eqs. (21), (22) and (24) depends 
on the factorization scheme convention adopted for the relation (22) between the NS 

^For the perturbative expansions of the QCD-/3-function and the anomalous dimensions we use the 
expansion parameter as /in (see, e.g., for detailed expressions). Note also that we have omitted the 
'A' for the (polarized) Jffsi^) since this quantity is trivially equal to its unpolarized counterpart H, [l9| . 
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structure function and the relevant NS quark densities beyond the leading order. Since 
the Qi^NS are physical, i.e. measurable, quantities and the 7Ars(^) is convention- 
independent, it becomes evident from eq. (21) that the scheme dependences of AC*(n) and 
A7j!^_5(n) cancel each other such that the combination 2AC*(?7,) + A7]^^(n)/2/3o is scheme 
independent flSl. Needless to say that removing all NLO quantities (AC* A7]^^, jSi) in 



eqs. (21-24), we recover the LO results of eqs. (9-14) with the quark density combinations 
evolving according to the LO (7^^^ |19|) NS evolution equation. 

Both the essential ingredients for the NLO calculation, ACg{n) and A'~fj^g{n), can 
be found in the literature. To facilitate the further discussion, let us first turn to the 
first moments, n = 1, which are of particular interest in the polarized case. As was 
discussed in ref. in the framework of the OPE, the operator corresponding to the 
first moments Ay43(l, Q"^) and Ay48(l, Q^) is nothing but the NS axial current which is a 
conserved quantity and thus has vanishing anomalous dimensions to all orders, which in 
particular means A7j^^(l) = 0. Furthermore, the value of the first moment of the Wilson 
coefficient ACl{n) for was found in |20|, |21| to be AC^i(l) = -3C^/2, giving 



—3Cf for the scheme independent combination 2AC*(1) -|- A7^5.(l)/2/3o and, according 
to eq. (22), leading to the factor (1 — Us/ti) in the NS sector of gi. Of course, both 
ACg(n) and A'-)\jg{n) depend on the regularization scheme adopted in their calculation^, 
and different schemes will in principle give different answers even for the first moments 
AClil), A7^s(l), though still respecting the condition 2ACl{l)+ A-i]^s{l) /2Pq = -SCp. 
However, the conservation of the NS axial current dictates the vanishing of A'~flfg{l), and 
hence the value AC^{1) = —SCp/^, which means that a scheme transformation has to be 
performed if these results are not automatically respected by the regularization scheme 
used. 

Let us briefiy list the results obtained for ACg(l) (to be calculated in the process 
7*g — * q{g) to 0{as)) using various regulators also used previously in the corresponding 
calculations in the unpolarized case. The result of ref. ||20|, AC^(l) = — 3Cir/2, was 
found using massless but off-shell incoming quarks and on-shell outgoing gluons. The 
same result for ACg{l) is obtained p2| if one uses massless on-shell quarks, but off-shell 



gluons {k > 0). Turning to dimensional regularization in MS, the result depends on how 
Needless to say that this has to be the same in the calculations of these quantities. 
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the Dirac matrix 75, appearing due to the projector on the quark's hehcity, is treated in 
n 7^ 4 dimensions. The prescription of a totally anticommuting 75 by Chanowitz et al. 
m yields |2|, |25| again]] ACl{l) = -3Cf/2. The same result is obtained [||] in the 
closely related 75 scheme of ref. |^^, taking the 7*g vertex as the 'reading point' to be 
defined in that scheme. However, when using the original scheme of 't Hooft and Veltman 
and Breitenlohner and Maison |]2D| (HVBM), or the equivalent prescription of refs. 



3T| , one obtains p5| , p2 



ACl{l) = -'^-Cp , (25) 

(naively) corresponding to a correction (1 — Tas/Sir) in the NS sector of gi and to a non- 
zero value for the anomalous dimension, A7^^(l), in contradiction to the conservation 
of the NS axial vector current. Finally, the same result, AC^(l) = —7Cf/2, is obtained 
for massive on-shell quarks {rriq 7^ 0) in the process 7*(f — * q{g). For completeness, 
we list all the results for ACg{n) for arbitrary Mellin-n in the Appendix. Comparing 
with the corresponding results [33-37] for the Wilson coefficient Cg{n) for the unpolarized 
structure function F2/2X in the various regularizations, one finds that all the AC^{n) with 
the property ACi(l) = -3C^/2 satisfy 

Cj(n) - AClin) =CJ- + (26) 

(corresponding to Cg(z) — ACg{z) = Cf(1 + z) in Bj0rken-a; space). This implied regu- 
larization scheme independence of Cg{n) — AC^{n) can be understood as follows: As a 
consequence of the factorization theorem, the difference C^^ {n) —2Cg{n), where C^^ {n) 
are the n-moments of the 0{as) (NS) quark corrections to the unpolarized Drell-Yan 
process qq l*{,g), has to be the same in any scheme (see, for example, 0)- The 
same is true (see [|^, ^) for the difference AC^^{n) — 2ACg(n) with the 0{as) quark 
corrections —AC^^{n) to the polarized Drell-Yan process qq l*{g)- On the other 
hand, the annihilating quark lines in this process trivially give AC^^(ri) = C^^ {n) if the 
regularization scheme used in the calculation of AC^^(?t,) respects chirality conservation. 
It then automatically follows that Cg{n) — ACg{n) is also the same in all such scheme^. 

"^The same result in dimensional regularization was found earlier in pl[ without specifying the 75- 
prescription. 

^Alternatively, one can see the expected scheme invariance of C^(n) — AC^{n) from the fact that 
p3| , Q C^{n) ^ C'q(") ~ Cf + ;7q:Y^ is scheme invariant (where C^{n) is the coefficient function of 
the unpolarized structure function F3) and from the similar appearance of F3 and gi in the hadronic 
tensor. 
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In contrast to eq. (26), we have for the calculation in the HVBM scheme and for the 
niq ^ calculation (which gave ACg(l) = —7(7^/2) 

CJ(„) - AC'(,,,) = (- + —) + I ^^^^ ^ ^ (27) 

According to our previous observations, these regulators then necessarily break the rela- 
tion ACf'^(ra) = C^^ {n\ i.e., break chirality. In fact, it was shown in ref. p5[ that 



1 

n n + 1 



in the HVBM scheme, showing how the terms ~ (1/n — l/(n + 1)) in eqs. (27), (28) 
cancel out in the difference AC^^(n) — 2ACq(r2), but individually break the relations 
ACg^^(n) = C^^{n) and Cl{n) - ACi(r2) = Cj.(l/n + l/(n + 1)). Furthermore, the 
terms ~ (l/n — l/(ra + 1)) in the HVBM scheme originate |2^ from a configuration where 
the incoming and the outgoing particles in the process 7*g qg become coUinear, and 
thus should rather be understood as part of the polarized (NLO) quark densities. Finally, 
as far as the massive calculation is concerned, the term ICpIn after the curly bracket in 
eq. (27) can be traced back to have its origin in a chirality breaking term ~ which 
survives the eventual limit — > since it happens to be multiplied by a double-pole 
term. Having found the origin of the additional terms in eq. (27) which lead to AC^(l) = 
— 7Cf/2, we expect that similar terms would be present in the I^'y^^in) when calculated 
in the HVBM or the massive scheme, such that scheme transformations, by means of 
2AC^(n) + /\'-)]qg{n) /2(3o = invariant, could be performed to eliminate these terms from 
both ACg and A'~flfg{n). Hereby one would obtain the correct values AC^(l) = —"iCp/'i 
and = (as dictated |20| by the the conservation of the NS axial current), and 



the relation C^iji) — AC^(n) = Cp{l/n + l/(n + 1)) would be restored in each case. 



To complete the discussion of the Wilson coefficients AC*, let us specify our final 
choices for the coefficients for gi, g^, and (74. Since, as we will see below, the anomalous 
dimension A'~fl^g{n) is known within dimensional regularization in the MS-scheme, we 
have to choose the Wilson coefficients accordingly. This means that the coefficient in 
(A. 4) ET], 123, pH] (or the one in (A. 5) after elimination of the chirality breaking term 



~ {1/n — l/(r;, + l))) is the relevant one for g[^'"^'\ It turns out [P2| , PB| that the coefficient 
is the same if electroweak contributions to gi (e.g. gi^) from transitions between massless 
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(A = 0) quarks q q' are considered. For the corrections to the structure functions g^, 
gi/2x one findsQ |2|, |8| 



ACl{n) 



^ \n n + 1 

^ n + 1 



(29) 



One notes the striking similarity to the relations between the quark corrections to the 
unpolarized structure functions F3, Fi, F2/2X which is readily explained by the similarity 
of the corresponding hadronic tensors. Equation (29) shows that the Callan-Gross-type 
relation = 2xg^ mentioned earlier is broken even for massless quarks beyond the LO. 
However, unlike its unpolarized analogue, Fl = F2 — 2xFi, which (in the singlet case) also 
receives contributions from g'/won-induced 0{as) corrections, g4 — 2xg3 = is broken only 
by quark-induced corrections (eq. (29)) even in the singlet case, since corrections from 
incoming gluons cancel out for massless produced quarks |]12 . 



In the case of the CC transition s — c we again have to take into account the mass 
of the charm quark which has an influence on the coefficient functions. For this purpose 
we have calculated the contribution of the process W~^s c{g) with mc 7^ to gi, g^, g^^ 
in MS dimensional regularization, following the techniques developed in The results 
of our calculation can be found in the Appendix. It should be noted that the expressions 
have a smooth limit ml/Q'^ ^ (A ^ 1), in which they reproduce eqs. (A. 4), (29). From 
eqs. (A.6-A.9) we immediately read off the 0{as) corrections to the functions fi{X,n) 
introduced in eqs. (11-14). For the first moment, n = 1, these functions then read in 
NLO: 



/i(A,l) 

/3(A, 1) 
/4(A, 1) 



Afl-3^C^^ln(l-A)^ 

\ 271 



A 

1-A 



A 
2 



-1-3A In(l-A) 



(30) 



The last result that /4(A, 1) receives no 0{as) corrections also holds for the corresponding 
function for the unpolarized structure function F2 |3D], where it is in accordance with the 
Adler sum-rule We emphasize that, similar to the unpolarized case 0, our 



^Eqs. (29) are actually independent of the regularization scheme chosen even in schemes where 
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results (A.6-A.9) for the contribution of the transition s — c to the spin-dependent 
structure functions would enable a determination of the proton's MS polarized strange 
quark distribution via a detection of charmed final states in polarized CC DIS. 

The (de facto) regularization scheme independence of the relation C^(n) — AC^{n) = 
CF{l/n + l/{n + l)) immediately implies that the scheme-dependent parts of the polarized 
and the unpolarized NS anomalous dimensions A'ylfg{n), 'ylfsi''^) equal each other in 
all schemes. Even more, as was first observed in |^0[ and recently established in more 
detail in the full expressions for A'yj^g{n) and 'jlisi^) exactly identical. This 
statement is correct in all regularization schemes, provided one has taken care to warrant 
ACq(l) = —3Cf/2 in the scheme used, eliminating, if present, chirality breaking terms 
as explained above. 

There is, however, another subtlety involved in the equality of A'yj^g^n) and 'ylfsi^)'- 
As is well-known [44-46], there is no analytical continuation of the unpolarized 7Ars(^) 
to arbitrary n, needed for the transformation from Mellin-n space into Bj0rken-x space, 
that reproduces the results for 7jV5(^) for all integer values of n. This is not unexpected 
since the OPE, first used to derive 'J^si''^) dimensional regularization [Q, gives 

only an answer for even n if the moments of the NS contribution to the e.m. structure 
function F2/X are considered, hereby artificially excluding odd values of n. Therefore, 
the analytic continuation of 'yj^fsi''^) only has to correctly reproduce the results for even 
values of n. On the other hand, as was shown in [|^|, odd values of n are relevant in 



the OPE for the NS combination F!^' /x — F!^"^ /x or for F^^ + F^^ , meaning that in this 
case the analytic continuation of 'y^si^) to reproduce the results at these values. 
These observations fit nicely and consistently to parton model considerations, where the 



NS quark combinations q — q' and qv = q ~ q can be easily seen to evolve pi| , |48|| with 
P+ = Pqq + Pqq and P = Pqq — Pqg, rcspectivcly, which are different beyond the LO, where 
Pqq, Pqq Bie the q ^ q and q q NLO NS splitting functions with flavor-non-diagonal 
contributions subtracted The explicit calculation of Pqq, Pqq shows that their 

Mellin-n moments satisfy[] 

7Nsin) = PM + i-irPM , (31) 



''For simplicity we have normalized the P± relative to 7^^^. 
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which means that the analytic continuation of 'yj^si''^) "which reproduces the values of 
iNsi''^) for evenn equals the combination P+{n) for arbitrary n, whereas the other analytic 
continuation of 'jlfsi''^)^ which is correct for odd n, corresponds to P^{n). In this way, the 
parton results of provide the rule for the analytic continuation of the OPE results. 
The essence of all this is that the moments of the combination Fg^^/x — Fg^^/x, or, more in 
general, the unpolarized y43(n, Q^), A^inyQ"^) (defined in analogy to eq. (7)) evolve with 
P-^-{n), whereas, e.g., the moments of /x — /x, F^^ + F^^ (which consist of pure 
valence, Uv{n,Q'^), dv{n,Q'^)), evolve with P^{n). 

The important difference in the polarized case is that the relevance of even and odd 
n in the OPE and for the analytic continuation is reversed here. As was shown in 0, ^ , 
odd n contribute in the OPE analysis to the combinations i^g^^^''^-'^ — g^^^'^'^){n,Q'^), 
{qT + Q^)) {q'z' ~^ 5'3")('^? Q^)) {qT ~ gT^){T^ ~ Q^)^ whereas even n are relevant, 

e.g., for {g7 - g7){n^Q^), {g? + 97){n,Q^), {g? + g7){n - l,Q^). In terms of the 
polarized NS quark distribution combinations this means that AA^{n^Q'^)^ /S.A^{n^Q'^) 
(as defined in (7)) evolve with P^{n) and the polarized valence densities Aut,(?T,, Q^), 
Ad^{n,Q^) with P+{n) This situation is summarized by the relations APgg = Pgg, 
APqq = —Pqq for the polarlzcd analogues, /S.P^q\ of Pq~q. 

4 Numerical results 

We are now equipped with all ingredients for a consistent NLO analysis of the spin- 
dependent NS structure functions. Let us consider the first moment of the NS combina- 
tions in eqs. (10,15-20). To begin with, we recall that the first moment -P-(l) = 7Ars'(l) = 
/S.'-j\jg{l) vanishes |^^. In the unpolarized case this is in accordance with the Adler sum- 
rule 1^ and the conservation of the number of valence quarks. In the polarized case it 
means that /S.A^{n,Q'^) and /S.A^{n^Q'^) for n = 1 correctly do not evolve with Q^, as 
required by the conservation of the NS axial current (see above): 

AA3,8(l,g') = AA3,8(l,Q^) . (32) 
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In contrast to this, the first moment of P+ is non-zero, which means that the first moment 
of the polarized valence densities evolves with beyond the LO: 

A,„(1,Q=) = + A,.(1.QS) , (33) 

with II P+(l) = ^CF{CF - Ca/2)(-13 + 12C(2) - 8C(3)) ^ 2.5576, and where we have 
used eq. (24) and the fact that the first moment of the LO NS anomalous dimension 
vanishes, 7Ars(l) = P, |19|. This yields the following sum-rules for the first moments of 



the polarized NS structure functions to 0{a^: 

^^ep(e.™.) _^e.(e.„..)) (^^ g2) = \ {x- ^ii^] A^ll, Qg) (34) 



^ ' A\ 37r 47r 2/?o / V3 3 / ^ ' ^° 

(35) 

2^/ (0, g ) = ^ (^1 + ^ ^ j (^3 Am. + -A4 j (1, Qo) (36) 

- .r) (1, Q^) = (l - ^ + (^^^ _ (3,) 

[97 - 97) (1, Q') = (l - ^^^) A^3(l, g^) (38) 
l{g7-9r)iO,Q')=^A,il,Ql) (39) 

(.f + 97) (1, Q^) - - (l - ^ + + A,,) (1, Ql) 

(40) 

^ (.f + 97) (0, g^) = - (l + ^^m^^^^^ (A.. + A4) (1, Ql) (41) 

2 + /^(A 1) ' ^1 ^ ^ 

etc. It should be noted that eq. (20) receives singlet contributions beyond the leading 
order, therefore we have not written down this equation any more. Eqs. (34-42) show how 
in principle measurements of the first moments of polarized NS structure functions even 
at large can serve to independently determine the combinations (Am. ± A(i.)(l, Qq), 
AA3{l,Ql) and AAs{l,Ql). This is particularly interesting considering the question 
raised earlier of which combination of polarized parton distributions can be related to 
the F, D values measured in baryonic /3-decays. To simplify the discussion, we follow 
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Q / ep(e.m.) en(e.m.)\ ^ ( n^'P \ n"' 

9[9i +9i j-o , .a .A9i +9i 



the recent NLO analysis to assume that at the low input scale Qq = 0.34 GeV^ 
(= filfLO [EH) can neglect any effects of SU{2) isospin breaking in relating /9-decay 
matrix elements of charged and neutral currents as well as SU (2) j breaking in the proton's 
polarized sea. We then have Ay43(l, Qg) = (Am^ — Ady){l, Ql) = F + D, and the rhs of 
eqs. (34,37-39) are completely specified, leading to unique predictions for the combinations 
of structure functions on the Ihs in NLO of QCD. The first of these is of course the well- 
known Bj0rken sum-rule |Q to The results for eqs. (34,37-39) are displayed 



in Fig. 1 as functions of Q^, where we have used ^Ms^^ ~ '^^^ MeV. To account for 



SU (3) breaking effects we parametrize the input quantities appearing on the rhs of eqs. 
(40-42) in the following way: 

AAs{l,Ql) = (3F-D)(l-ei) (43) 
(Am, + A4)(1,Q2) = ^^(3ir_/)) , (44) 

1 — 62 

which yields {3F - D) /2{l-e2) + {F + D)/6 for the combination (2Am^/3 + A4/3)(l, Qg) 
in eqs. (35,36). Eqs. (43,44) are general enough to take into account all possible sources 
of SU{3) breaking: ei determines the deviation of the first moment AAs{1,Qq) from the 
value 3-F — D obtained from hyperon /3-decays. Such a deviation will occur if the use of 
SU{3) symmetry for relating the matrix elements of charged and neutral axial currents 
is not justified. In this case, ei could be significantly different from zero, even such that 
only the valence quarks contribute to 3-F — -D • This possibility is taken into account 
by the parameter €2 which would vanish in the latter case. From the definition of AAg 
one furthermore sees that ei and €2 together determine the amount of SU{3)f breaking 
in the proton's polarized sea: 

Auy + Ady 

Fig. 1 shows our predictions for the NS structure functions of eqs. (35,36,40-42) for the 
conceivable choices [^] ei = 0, £2 = 0.105 and ei = 0.40, €2 = 0. It becomes obvious that 
the effects of changes in the are larger than the present experimental 4%-uncertainty ^ 
in the value for 3F — D (see eq. (4)) and that therefore a measurement of the quantities 
shown would help to decide about the amount of SU{3) breaking. In particular, the 
parameter €2 could be determined in NC,CC experiments with polarized beams at HERA 



*Note that actually the corrections to 0{al), 0{al) to this sum-rule are known pi| 
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g via a measurement of f^f '"^^(0, Q^) or ((^f + g'^^){0, Q^) = {gl ''^^ + gf^'^^M Q') 
(or their (73- analogues). Using the full Mellin-n-dependent expression for ACg{n) from 
eqs. (29), (A. 4) in eq. (22), we can obtain NLO predictions for the Bj0rken-a; dependence 
of the latter structure functions: 

<7f'^^(n - 1, Q') ] / \ r i (|A«. + |A4) in, Q^) 

(^r•''^ + ^f^''^'')(r^-l,g^)J \ 27r ^""^ I ("2) (A«. + Aci.) (n, Q^) , 

(46) 

where, again, the polarized valence quark densities are to be evolved according to eq. 
(24) with the correct analytic continuation P+{n) of the 'j'^si''^) found in |5( 



The results for gT'^^ix,Q^) and {gl ^'^^ + gl P'^^)(x,Q2) at = 1000 GeV^, found 
after Mellin-inverting eq. (46), are shown in Fig. 2, where for the polarized input valence 
densities Aqy{x,Ql) at Ql = 0.34 GeV^ we have used the two sets determined in ref. 
P9|| . Both sets give a very good description of all existing data on deep- inelastic spin 
asymmetries in the valence region x ~ 0.2, but they differ in the assumptions made about 
the role of SU (3) / symmetry breaking effects and therefore have different first moments 
||49|| , corresponding to the ei, €2 values used in Fig. 1. Thus the variation in the results 
shown in Fig. 2 for the different sets of polarized valence input densities reflects the 
present theoretical uncertainty in the predictions. Conversely, Fig. 2 shows that also a 
measurement of g'^'^'^{x,Q'^) and ((74 ^''"^ + gl ^''"'"){x,Q^) for x < 0.2 at HERA could 
help to shed light on the importance of SU (3) / symmetry breaking. 

The different evolution of the polarized valence quark densities and the combination 
AA3 beyond the LO induces a dynamical breaking of the SU{2) f symmetry in the proton's 
polarized sea^. Eq. (24) and our considerations concerning the analytic continuation of 
ihsi^) predict 

(«.(Q')-«.(Qg)) (P-(^)-nH) 



2 (A?i - Arf) [n, 



47r 2Po 

2\\ 7^sW/2/3o 



xf^4^) {Au^~Ad,){n,Ql) , (47) 

where we have again assumed the nonperturbative input at Qq to be SU{2) /-symmetric, 
Au{n, Ql) = Ad{n, Qq). Using again the polarized input valence distributions of ref. 



^There is also a dynamical breaking of SU{3)f symmetry in the sea induced by AA^. We do not 
pursue this effect since it is most probably dominated by the SU (3) f breaking in the nonperturbative 
input (see eq. (41)) due to the larger strange mass. 
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at Qq = 0.34 GeV^, we obtain a prediction for {Au — Ad)(x, = lOGeV^) which is shown 
in Fig. 3. Both sets of polarized input valence densities considered in |^ lead to entirely 
indistinguishable results, since only the input combination {Au^ — Adv){x,Ql) is needed 
here whose first moment is fixed by the value F + D (eq. (3)) in both cases. It should be 
noted that such a dynamical breaking of SU (2) f symmetry in the sea induced by two-loop 
evolution was considered in the unpolarized case in ref. EH], where it was found to be very 



small. The results in the polarized case differ in sign (due to the interchange P+ ^ P_) 
and slightly in magnitude (due to the polarized valence input instead of the unpolarized 
one) from the unpolarized results. However, the relative effect of SU{2) f breaking is much 
larger in the polarized case since the polarized sea densities are probably much smaller 
than their unpolarized counterparts. Even so, when taking the first moment, one obtain^ 

2(A^-Acl)(l,g^ = 10GeV^) = _ ("^(Q^) " "-(^o)) ^jF + D) 

^ ^ 4:71 

^ 0.006 , (48) 

which means that unless the input at Ql strongly breaks SU{2)f symmetry the effect 
of the breaking is probably small compared to the size of |A'u(l, Q^)|, |Ac/(l,(5^)| which 



might well be of the order ~ 0.05 |^9|. It is straightforward to introduce parameters ^i, 82 
in analogy to ei, 62, which would parametrize genuine SU{2) breaking effects in the first 
moment of the polarized sea and in the relation between charged and neutral axial current 
/9-decay matrix elements. Measurements of the first moment of the structure functions in 
eqs. (34,37-39) (see also Fig. 1) would then allow to determine these parameters and to 
pin down SU (2) breaking effects. 



5 Summary and conclusions 

We have performed a detailed study of spin- dependent non-singlet structure functions 
in the framework of a complete and consistent NLO QCD calculation. Our analysis is 
based on a careful discussion of the calculation of the 0{as) corrections to the structure 
functions, in which we have examined the regularization scheme dependence of the NS 
coefficient function AC^ for gi with respect to the constraints imposed by axial current 

-"^"This number depends quite crucially on the value chosen for the input scale Qq. Taking, e.g., Qg = 1 
GeV^ the result in eq. (48) is reduced by a factor 3. 
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conservation. We have also shown how to correctly take into account the two-loop evo- 
lution of polarized NS quark combinations. A further ingredient of our study is the full 
inclusion of the charm mass effects in the charged current s — > c contributions to polarized 
electroweak structure functions. 

Our numerical analysis has revealed that conceivable measurements of spin-dependent 
NS structure functions at HERA or in scattering experiments off polarized nucleon 
targets would serve to improve our understanding of the relations between the first mo- 
ment of NS combinations of polarized quark densities and the F, D values extracted from 
hyperon-/5-decays, and would also shed light on SU{2)f, SU{3)f breaking in the nucleon's 
polarized sea. Finally, we have also shown that the latter symmetries are dynamically 
broken by NLO evolution in the NS sector. 
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Appendix 



In this Appendix we list the results for the polarized coefficient function AC^{n) using 
various regulators in its calculation from the process 7*g —>■ q{g)- In all cases we have 
chosen to just subtract the collinear pole contribution, which is then factorized into the 
(bare) quark distributions. The singular terms are of the forms 'y%g{n)ln{Q'^/\m'^\) (if 



some mass or off-shellness y\m? \ is used as the regulator) or 'j%g{n){—l/e) (in dimensional 
regularization in the MS-scheme). Since the first moment of 7Ar5(^) vanishes [P, 19|, the 
pole contribution drops out from the more important first moment anyway. The results 
in Mellin-n space below can be easily transformed into Bj0rken-x space with the help of 
the Appendix of ref. p4|. 



Off-shell massless quarks, on-shell gluons: This calculation corresponds to the one 
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of ref. PDI, but our results slightly differ by a term Cf{—^ + 2/{n + 1)) (which 
vanishes for n = 1) due to the specific operator normalization chosen in ||2^ (see 
the Appendix A. 2 of ref. pO[ for details). For the result in the unpolarized case {F2) 

2 



see 



2 2 

2n ' n + 1 



3 



o n 1 '^1 



(A.l) 



On-shell massless quarks, off-shell gluons: The AC^{n) for this calculation can be 



obtained from |22|. For the unpolarized case see |35|. 



ACi(n) 



9 



4 2n n+1 it? (n + l)^ 
1 



n 



(A.2) 



On-shell massive quarks, off-shell gluons: In this regularization we obtain: 

2 1 2 



ACl{n) 



5 5 

~2 ~ 2^ ^ n + 1 ^ ^ ~ (n + 1)2 

7 1 \ A 1 A 1 A 1 1 

o + —f 

2 nin + 



Note that ACg(l) = —lCpl2 in this scheme. See refs. 
unpolarized result. 



(A.3) 



for the corresponding 



Dimensional regularization: Using the 75-prescription of ||23| (or its more systematic 



and consistent generalization ||2^) one obtains in the MS-scheme |2 



AC\[n) 



9 1 1 
'2^2^^ n + l 



1 



)n 1 "1 1 -J 1 

y--2y- + 2y-y- 



(A.4) 



The same result in dimensional regularization was found earlier in |^ without 
specifying the 75-prescription. However, using the original scheme of 't Hooft and 
Veltman and Breitenlohner and Maison (or the equivalent one of refs. 
[0, |3ll]), one finds ^ an additional term 

1 1 



(A.5) 



which leads to AC^(l) = —7Cp/2. For the unpolarized case see ^ 
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We finally present our results for the coefficient functions AC* for gi, g^, gi/^x for 
the transition s — > c, fully taking into account the effects due to the charm quark mass. 
The calculation was performed in MS dimensional regularization in the 75-scheme of ref. 
27|, choosing the axial vertex as the reading point. Our Bj0rken-a; space results for AC*^, 
ACg, ACg fully agree with those of ref. for the unpolarized h^^q hi^q, h2,q (for F3, 
Fi, F2/2X), respectively, after eliminating an error in the coefficient A2 in that paper 
which should read Ka instead of Ka/^. The differences AC^ — AC^, AC^ — AC^ which 
are regularization scheme independent, are in agreement with the results of for the 
corresponding differences in the unpolarized case. We note that the results of ref. [H^ 
seem in slight disagreement with both [p9| (even after correction of the above mentioned 



error) and |^ and also with our calculation in this respect. Here we present the Mellin-n 
moments of our results. For this purpose it is convenient to present the moments for the 
differences AC*(n, A) — AC*(n), where the ACg{n) are the (usual) massless coefficient 
functions given in eqs. (29), (A. 4), and A = Q^/{Q^ + ^c)- Defining the sum 



and 39 



Sa{n, A) = A 



Ka{X) = 



ln(l - A) + E 



1-A 
A 



ln(l - A) 



we find: 



ACl{n,\)-ACl{n)=C, 



1)1 



1=1 



in — \)in + 2) , „ , , , ^ , , , , nin 
+ ^ . . \ ' (5o(n, A) - 5i(n, A)) ^ 



A 




2{n + l) A 



Si{n,\) 



--2y 

1 



In A 



ACj-ACn(n,A) 



ACl 



[AC, 



AC: 



[n. 



A) - (AC, 



ACl 



AC: 



in 



Cf 
Cf 



1-A 1 



A 



\\2 

2-^5i(n,A) 



A n + 1 

(1-A^)(1-2A) 



(A.6) 
(A.7) 



Ka{\) - 



A2 



Si{n,X) 
(A.8) 



n X{n + 1) y 

The last term in eq. (A.6) which contains the LO 7Ars(^) 
chooses the scale as the factorization scale It should be noted that, like in the LO 



T9|] is introduced if one 
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(see eqs. (11-14)), an additional factor of A" (A" ^) is needed to calculate the contribution 
to the structure functions gi{n, Q^), g^in, Q^) {g4,{n — 1, Q'^)/2). 
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Figure Captions 



Fig. 1 Predictions for the Q^-evolution of tlie first moments of the various NS combina- 
tions of polarized structure functions as given in eqs. (34)-(42) for two conceivable 
choices of SU{3)f breaking parameters ei, €2 in eqs. (43), (44). The input scale for 
the evolution, Ql = 0.34 GeV^ was chosen according to ref. [^, and ^^(Q^) was 
calculated from eq. (23) with from |^ . 



Fig. 2 Predictions for the NC and CC non-singlet structure functions (cf. eq. (46)) 
(74^'^*" (x, Q^) and gT''^'^{x,Q'^) = (gl + (74 (x, Q^), respectively, as mea- 



surable in a future polarized e~p/e^p collider mode of HERA |jT6|. For the predic- 
tions we have used the two sets of polarized input valence densities suggested in 
which correspond to the SU{3)f breaking parameters introduced in Fig. 1. 

Fig. 3 Prediction for the dynamical SU{2)f breaking of the proton's polarized sea 
(^Au — A(tj {x , Q'^) at = 10 GeV^ according to eq. (47). The nonperturba- 
tive valence input (Am^, — Ad^) (x, Ql) at Ql = 0.34 GeV^ was taken from the 
analysis in ref. For comparison the dashed line shows the averaged sea den- 

sity —Aq{x,Q'^) = — (^Au{x,Q'^) + Ad{x,Q'^)^ /2 determined within the 'standard 
scenario' of ref. ET 
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T — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 

eq.(39) 
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